Andreev interferometer with three superconducting electrodes 



Artem V. Galaktionov 

I.E. Tamm Department of Theoretical Physics, P.N. Lebedev Physical Institute, 119991 Moscow, Russia and 
Laboratory of Cryogenic Nanoelectronics, Nizhny Novgorod State Technical University, 603950 Nizhny Novgorod, Russia 

Andrei D. Zaikin 

Institut filr Nanotechnologie, Karlsruher Institut fiir Technologic (KIT), 76021 Karlsruhe, Germany and 
I.E. Tamm Department of Theoretical Physics, P.N. Lebedev Physical Institute, 119991 Moscow, Russia 

Leonid S. Kuzmin 
Chalmers University of Technology, Gothenburg, Sweden and 
Laboratory of Cryogenic Nanoelectronics, Nizhny Novgorod State Technical University, 603950 Nizhny Novgorod, Russia 

We develop a quasiclassical theory of Andreev interferometers witii three superconducting elec- 
trodes. Provided tunneling interface resistance between one superconducting electrode and the 
normal metal strongly exceeds two others, significant current sensitivity to the external magnetic 
flux is observed only at subgap voltages. If all barrier conductances are comparable, multiple An- 
dreev reflection comes into play and substantial current modulation can be achieved in both subgap 
and overgap voltage regimes. Our analysis reveals a large variety of interesting features which can 
be used for performance optimization of Andreev interferometers. 



I. INTRODUCTION 

Andreev interferometers are often regarded as poten- 
tial rivals to Superconducting Quantum Interference De- 
vices (SQUIDs) for a number of applications with the 
possibility to achieve higher sensitivity and read-out 
speed. Such kind of applications range from studying of 
switching dynamics of individual magnetic nanoparticles 
to read-out of superconducting qubits. 

Almost two decades ago Petrashov and co-workers 
performed experimental analysis of magnetoresistance 
of mesoscopic normal-superconducting (NS) hybrid 
structures"'^. This analysis revealed conductance modu- 
lation greatly exceeding universal conductance fluctua- 
tions. It was perceived, that this modulation should be 
attributed to the effect of the phase difference x between 
superconducting elements of the hybrid structure influ- 
encing the process of Andreev reflection (see Refs. 
for a review). Since the normal bars in the cross-like 
diffusive structures studied in experimentsi exceeded the 
superconducting coherence length the characteristic 
energy scale for magnetoresistance modulation should be 
set by the Thouless energy et/i — D/X'^, where D is the 
diffusion coefficient and X is the bar length in the cross. 

In a diffusive NS structure with good transmission of 
metallic interfaces the magnitude of the magnetoresis- 
tance modulation is expected not to exceed few percents. 
This estimate is directly related to the so-called reen- 
trance effect in the corresponding NS structures reaching 
the maximum value below 10 percent^i^. If, however, 
tunnel barriers are present in the NS system its conduc- 
tance can change by much higher values^ . This is because 
Andreev conductance of a tunnel barrier at the NS inter- 
face is much smaller than its normal state conductance 
while a diffusive connector has the same low tempera- 
ture Andreev conductance as its normal one. Hence, one 



could expect that magnetoresistance modulation could 
also be much more pronounced in hybrid NS structures 
which contain tunnel barriers. This feature was indeed 
demonstrated experimentally by Pothier et al^ who ob- 
served the maximum-to-minimum resistance ratio as a 
function of x to be as high as ~ 5. Physical insight into 
the phase dependence of the Andreev conductance in hy- 
brid NS structures is provided by generalization of the 
Kirchhoff rules worked out by Nazarov'''. 

In the experiments discussed so far at least one of the 
external electrodes was in the normal state. Yet another 
option to fabricate the interference device is to keep all 
available electrodes superconducting. In this way one 
would be able to reduce dissipation. Current harmonics 
of the current would be generated in this case, which are 
multiples of the Josephson frequency. However, they can 
be filtered out and the average current can be measured 
which should reveal a dependence on the phase differ- 
ence X- Experiments with such structures were recently 
performed by Meschke et al^ , and the corresponding the- 
oretical analysis was developed in Ref. |9|. 

In this paper we study theoretically a different setup 
which is schematically depicted in Fig. [TJ This device 
consists of a disordered normal insertion embedded in- 
between three superconducting electrodes. Typical size 
of this normal insertion L does not exceed the supercon- 
ducting coherence length = \/D/ A, where A is the 
superconducting gap of the electrodes, D = vpl/'i and I 
is the elastic electron mean free path. At the same time, 
the normal metal size obeys the condition L ^ I. There 
exists a superconducting phase difference x between the 
electrodes 2 and 3 which is controlled by the external 
magnetic flux $ piercing the superconducting loop, i.e. 
X = 27r$/<l'oj where $o is the superconducting flux quan- 
tum. 

In what follows we will generally assume that inter- 
faces between normal metallic dot and superconducting 
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FIG. 1: Scheme of the setup: disordered normal insertion 
("dot") between three superconducting electrodes. The phase 
difference X2 — Xs — X is caused by the magnetic flux piercing 
the loop. 



electrodes are weakly transniitting and their normal state 
conductances Gi, G2 and Ga are supposed to be smaller 
than the dot conductance ~ ao-A/L. Here ao — 2e^DNo 
is the Drude conductivity of the normal metal, A is the 
typical contact area between the normal metal and the 
electrode and iVo is the density of states at the Fermi- 
surface per spin direction. This relationship between 
conductances assures the voltage drops only across NS 
interfaces, while there exist no significant voltage varia- 
tions inside the normal dot. Under these conditions our 
results will not depend on the particular shape of the nor- 
mal metal insertion. For example, one can equally apply 
our analysis to the setup displayed in Fig. [2] We will also 
ignore charging effects which amounts to assuming that 
all relevant charging energies remain much smaller than 
the corresponding Josephson coupling energies^ and, in 
addition, that tunneling conductances Gi,2,3 strongly ex- 
ceed the quantum conductance unit Gq — e^//*ii. The 
latter condition also allows to improve noise characteris- 
tics of the systems under consideration. 

In the analysis^i^ it was assumed the contacts between 
electrodes 2, 3 and the normal insertion are ideal i.e. 
highly transmitting) while the first electrode is connected 
to the dot via a tunnel barrier. Obviously in this case 
the conductances obey the condition G2 « G3 ^ Gi. A 
proximity-induced x-dependcnt minigap inside the nor- 
mal dot develops in this case. As a result, the current- 
voltage characteristics of such Andreev interferometer 
should be sensitive to the phase difference x which, in 
turn, can be controlled by external magnetic fiux $. 
Qualitatively the same features hold in the situation 
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FIG. 2: Another possible realization of interferometer with 
the Y-shaped normal metal insertion. 

studied here as well. On the other hand, there also ap- 
pear important differences. For instance, we will demon- 
strate that the properties of the Andreev interferometer 
essentially depend on the parameter 7 ~ tjjA, where td 
is the electron dwell time in the normal dot. Similarly to 
the case of diffusive SNS junctionsi^ this parameter ef- 
fectively controls the strength of electron-hole dephasing 
in our device. Another non-trivial feature of our struc- 
ture is related to the effect of multiple Andreev reflec- 
tions (MAR)i^ which becomes particularly important in 
the case Gi ~ G2 ^ G3 providing significant modifica- 
tions in the x-dependent current both in high and low 
voltage limits. The effect of MAR on the properties of 
Y-shaped Andreev interferometers was studied in Ref. [13 
for a model system with normal insertion substituted by 
a one channel quantum wire. Here we will address a real- 
istic configuration which includes many-channel diffusive 
conductors of an arbitrary shape. 

The structure of our paper is as follows. In section II 
we outline our general formalism which allows to evaluate 
quasiclassical electron propagators for the problem under 
consideration. Section HI is devoted to the analysis of the 
density of states in the normal metallic dot as well as the 
flux-dependent current across our structure in the regime 
Gi ^ G2,3. Non-perturbative MAR regime is studied 
in section IV. Our main conclusions are summarized in 
section V. 



II. QUASICLASSICAL ANALYSIS 

Our theoretical analysis is performed with the aid of 
Usadel equations (see, e.g., Ref. H) 

-Ddni9odng) + f^^ + ^f^+ (1) 

(-iA(i) + ieip{t)) g ~g {~iA{t') + ieip{t')) = 0, 

where g represent the quasiclassical propagators that are 
4x4 matrices depending on one spatial and two time 
variables 

9iR^i^n-[ g^iR,t,t') )■ 
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Here g^^^'^ stand for retarded, advanced and Keldysh 
components respectively, each of them forming a 2 x 2 
matrix. The electric potential is denoted by (p{R, t) and 
e is the electron charge. The product of propagators in 
Eq. ([!]) implies time convolution 

oo 

{giog2){t,t')^ [ dhg{t,h)g2{h,t')dti. 



The remaining matrices in Eq. ([T]) are defined as 



f, 



1 



Of,/''" V -1 / ' " V A y ' 



A 



where 



A 



A(H,t) 
-A*{R,t) 



(3) 



and A is the superconducting order parameter. Qua- 
siclassical propagators ([2]) also obey the normalization 
condition 



(gog) {hM)^5{h-t2). 



(4) 



In the immediate vicinity of a tunnel barrier Eqs. ([T]) are 
not applicable and the quasiclassical propagators evalu- 
ated on both sides of this barrier should be matched by 
appropriate boundary conditions^ 



2(7 ongdRg = ctt [9~,9+\ ■ 



(5) 



Here the subscripts ± label the propagators on the right 
and the left sides of the interface, n is the unit vector 
perpendicular to the interface (directed from side "-" to 
side " +" ) . The combination in the left-hand side of Eq. 
([S]) is continuous at the interface, the commutator in the 
right-hand side is denoted by brackets and ctt stands for 
conductivity per unit square of the tunnel barrier. 

It is convenient to Fourier transform the quasiclassical 
propagator with respect to the time difference 



g{R, e,t)= j dt' exp(iei')ff (-^' ^ + |' ^ " | 



(6) 



Integrating Eq. ([T]) over the normal metal volume V and 
employing the Gauss theorem one finds 

V [ef,,t] - I j dVDdk {gdk'g) = j dSk {gdk'g) ■ (7) 

With the aid of the boundary conditions ([5]) from Eq. ([7]) 
we obtain 

[Z,f\=^, Z = ^angn~ il^^Tz, (8) 

where the unperturbed quasiclassical propagator of the 
n-the electrode is denoted by ^„ and we also introduced 
the notation 



Gi + G2 + G3 ' 



n = 1,2,3. 



(9) 



Finally, we introduce the parameter 



7 ^ 



2az,V|A| 



D{Gi +G2+ G3) 



rD\A\ 



(10) 



which is a direct generalization of the analogous parame- 
ter defined for diffusive SNS junctions^^. This parameter 
effectively controls the strength of electron-hole dephas- 
ing in our system. 

Combining Eq. ([5]) with the normalization condition 
(|4|) we arrive at the expression for the quasiclassical prop- 
agator in the normal dot. It reads 



g = 



/Z2 



(11) 



We note that such an expression was previously discussed 
in Ref. [l^ in the context of the full-counting statistics 
and in Ref. [l^ in the context of MAR in SNS junctions. 



III. DENSITY OF STATES AND TUNNELING 
CURRENT 

Let us first consider the limit Gi <C G2,G3. In this 
case Eq. (|5]) yields 



Z = 



G2g2 + G3g3 
G2 + G3 



(12) 



The retarded and advanced components of the equilib- 
rium quasiclassical propagators of the electrodes 2 and 3 
are 



R,A 



e^^ = ±^{e±t6)^-\A\^ (13) 



while the Keldysh component is defined as g^ = g^F — 
Fg'^, where F{e) = tanh(e/2T) is the Fourier-transform 
of the function F{t) = —iT/ sinh[7rTt]. We also assume 
that there exists the superconducting phase difference be- 
tween the order parameters in the electrodes 2 and 3, i.e. 
in these two electrodes we define A = \A\e^'^^^. As we 
already pointed out, this phase difference is proportional 
to external magnetic flux $ piercing the superconducting 
ring. 

Considering the diagonal component of the matrix 
{g^ — g^)/2 and employing Eq. (fTTj) we recover the den- 
sity of states in the normal dot equal to 



n(e) — Re 



\/e^ - A2 ' 



where we define 



A = 



Here the parameter 



l + 7Vl-eV|A|= 



- lAL 1- 



4G2G 



(G2 



• 2 X 

G,y 2 ' 



(14) 



(15) 



(16) 
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FIG. 3: (Color online) The density of states defined in Eqs. 
(O-Hni with eg = 0.5|A|. The value of 7 is 1.0 for the upper 
plot and 10 for the lower plot. 
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FIG. 4: (Color online) The current-voltage characteristics 
evaluated from Eq. ([T8| for T = 0.1 A and 7 = 1. The 



this barrier reads 



2eR 



N 



den{e)nQ{e + eV) x 



(18) 



, , e+eV f 

X tann tann — , , 

2T 2T ' 



represents the minigap in the density of states of the 
normal dot provided the dephasing parameter tends to 
zero 7 0. In the symmetric case G2 — G3 Eqs. p^ - 
(fT6)) coincide with the corresponding expressionsii. 

In order to interpret the resuh (jl6p in a general case 
G2 7^ G3 we recollect that the maximum transmission 
value Tniax for the system of two tunnel barriers with 
transmissions T2 and T3 (both much smaller than one) is 
defined by the well known formula 

It follows immediately that Eq. (IT5|) just defines the po- 
sition of the Andreev level with energy 

^a{x) = |A|Y/l-T^axSin2(x/2). 

For non-zero values 7 the minigap Ag in the normal dot 
is obtained from the solution of the following equation 

A, = , (17) 

1 + 7.^1-A2/|AP 

i.e. in the limit of strong dephasing 7 3> 1 the minigap 
gets reduced as Ag w ^g/l- The corresponding density 
of states is exemplified in Fig. [31 

Let us now evaluate the dissipative current flowing 
across our device at a given voltage bias V . In the limit 
Gi <^ 6*2,3 considered in this section the voltage drop 
concentrates at the tunnel barrier between the first elec- 
trode and the normal metallic dot. The current across 



where Rn ~ 1/Gi and 

\e\e{\e\-\A\) 

no(e) — , - (19 

^ ' ^£2- |A|2 ^ ' 

is the standard BCS density of states in the first elec- 
trode, 9{x) denotes the Heaviside step function. Exam- 
ples of the current- voltage characteristics evaluated from 
Eq. p8)) are displayed in Fig. |4l We observe that pecu- 
liarities (spikes) on the I-V curve occur at voltages equal 
to eV = A - Ag, eV ^ A + Ag and eV = 2A (less pro- 
nounced). Here and afterwards A stands for the modulus 
of the superconducting order parameter. 

Typical curves displaying the phase dependence of the 
tunneling current (ITS)) at voltages eV > A are shown in 
Fig. [5] We observe that provided temperature is low 
and eV < 2A there exists a pronounced jump in the 
current which occurs as the phase x reaches the value 
corresponding to eV = A -t- Ag(x). For eV > 2 A the 
phase dependence of the current is monotonous, the cur- 
rent decreases with growing x implying " positive magne- 
toresistance" of our structure. The current modulation 
amplitude decreases rapidly as the voltage increases and 
becomes almost negligible already at eV ^ 3A . 

As long as temperature is not too low and the value 
exp(— A/2T) is not vanishingly small the X"dependent 
current is also observed at voltages eV < A. The corre- 
sponding plot is shown in Fig. [6] One observes a pro- 
nounced current peak which occurs at eV = A — Ag (x) . 

In order to complete this part of our analysis we should 
add that in order to avoid hysteretic phenomena in our 
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FIG. 5: (Color online) Phase-dependent tunneling current in 
a symmetric device at T = 0.1 A and 7 = 0.5. 




FIG. 6: (Color online) Phase-dependent tunneling current in 
a symmetric device at T = 0.1 A, 7 = 0.5 and eV = 0.7A. 



Andreev interferometer it is necessary to obey the con- 
dition 



£7£ ^ £0 

c 2t:'' 



(20) 



which is exactly analogous to that well known for stan- 
dard SQUIDs, see, e.g., Ref. JJ). Here £ is the loop in- 
ductance and Ic is the critical Josephson current in our 
device. The value Iq can roughly be estimated with the 
aid of the Ambegaokar-Baratoff formula for the Joseph- 
son current 



Is (x) = TT-^ sm X tanh — 



(21) 



with i?23 = + ^3^^- Strongest deviations from this 
formula occur for symmetric structures with G2 = G3 
provided 7 — >■ 0. In order to estimate these deviations in 
this case we can use the following transmission distribu- 
tion of conducting channelsi^ 



P(T„) cx 



1 



(22) 



Employing this distribution and evaluating the Joseph- 
son current, in the limit T we obtain 



A 

ei?23 



smx 



K sm — 
V 2 



(23) 



Here K{x) = £'^d(j)/^/l 



X sin (j) is the complete el- 
liptic integral. The maximum current value provided by 
Eq. (|23p is higher than that determined from Eq. (PT|) 
by a factor « 1.22, i.e. maximum deviations from Eq. 
(|2T|) are in the range of 20 percents. 



IV. BEYOND TUNNELING LIMIT: 
REGIME 



MAR 



In the previous section we restricted the analysis of I-V 
curves to the tunneling limit, i.e. evaluated the current in 
the lowest order in the conductance of the tunnel barrier 
with lowest transmission Gi <^ G2,3- If this condition is 
violated, e.g., all three conductances have approximately 
the same value Gi ~ G2 ~ G3, it is necessary to go be- 
yond the lowest order perturbation theory and account 
for higher order electron tunneling processes between su- 
perconducting terminals. In this case an important role 
is played by the mechanism of multiple Andreev reflec- 
tion (MAR): quasiparticles with energies below the su- 
perconducting gap A propagating inside the normal dot 
and suffering Andreev reflections at different NS inter- 
faces (i.e. being converted from electrons to holes and 
back) are accelerated by the electric field and eventu- 
ally leave the dot area as soon as their energies exceed 
A. Recently it was demonstrated theoretically that this 
process essentially influences the I-V curves of diffusive 
SNS junctions^^. The system under consideration here is 
more complicated since it contains three superconducting 
terminals. 

Let us set the electrostatic potential of the first elec- 
trode equal to —V. Then the quasiclassical electron prop- 
agator in this electrode reads 



gi(i,t')=e'^^^^*geq(t-t')e 



(24) 



where the equilibrium propagator geq{t) is defined in Eq. 
P^ . The general expression for the current flowing into 
this electrode through the tunneling barrier with conduc- 
tance Gi has the form 

/(i) = ^GiTr(f4t,5i]'')(i,i), (25) 

where the superscript K denotes the Keldysh component 
of the commutator. The electron propagator inside the 
normal dot g is given by Eq. , where the matrix Z is 
deflned in Eq. ([5]). One can also cast the expression for 
g to the form 

00 

t==;^y" dXK{\), K{X) = {Z + (26) 
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Provided all electrodes are in the normal state the above 
equations just yield Kirchhoff rules, i.e. 



Gi + G2 + G; 



■{G2{V2-Vi) + G3iV3~Vi)). (27) 



In what follows we will make use of this expression in or- 
der to normalize the corresponding results derived below 
for the superconducting case. 

It is convenient to rewrite the dependence of the qua- 
siclassical propagators (H]) on e and t with the aid of a 
series in multiples of the Josephson frequency, i.e. 



-2imeVt 



(28) 



m— — 00 



In this representation the matrix Z + iX has only m = 
— 1,0,1 components which will be denoted as i?_i(e), 
Ho{€,X) and Hi{e) respectively. The condition 



{Z + iX) o K ^ S(t-t') 



reduces to 



(29) 



Hoie + meV, X)Kxie,m) + 
Hi{e + (to - l)eV)Kx{e - eV, m - 1) + 
H-iie + (to + l)eV)Kx{e + eV, m + 1) = (5,„,o- 



Introducing K„i{e, X) = K\{e-\-meV, m), one can rewrite 
Eq. ^ in the form 



Hoie + 2meV,X)Kmie,X) + 
H,{e+{2m-l)eV)K,n-i{e,X) + 
H^i{e+{2m+l)eV)K^+i{e,X) = 



(30) 



In order to resolve this equation one can employ the 
ansatai^ 

K,n = S„iS,n-i ■ ■ . SiKo; m > 0, (31) 

Km — PmPrn+1 ■ ■ ■ P-iKq', TO < 0, 

with the aid of which Eq. (|30p yields recurrences relating 
Sm and Sm+i 

5„(e,A) = - [H^i{e + (2m + l)eV)Sm+i{e,X)+ 

Ho{e + 2meV, X)y^ Hi{e + {2m - l)eV). (32) 

Similarly, for to < 0, we obtain the following relationship 
between and Pm~i 

Pm{e, A) - - [Hi{e + (2m - l)eF)P,„_i(e, A) + 
Ha{e + 2meV, A)] H-i{e + {2m + l)eV). (33) 

At TO = we get from Eq. 

i^o(e,A) = [Hi{e-eV)P-i{e,X)+ (34) 
H-i{e + eV)Si{e,X)+Ho{e,X)y\ 
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FIG. 7: (Color online) Current-voltage characteristics at 02 = 
Qfs = 0.45, T = O.IA and 7 = 1. The values of x are 0, 1.3, 2 
and TV. The lower x value corresponds to the upper curve for 
eV > 2A. 




FIG. 8: (Color online) Phase-dependent current at 02 = as 
0.45, r = O.IA, 7 = 1 and = 1.3A. 



Finally, we impose the "boundary conditions" 

0, lim P„ = 0. (35) 

1— V— 00 



lim S„ 



Thus, proceeding numerically we set Sm, Pm = for some 
large \m\ and then employ the recurrences ([5^ and ([55)1 
in order to find 5*1, P_i and Kq. 

Averaging of Eq. (|25|) results in the following expres- 
sion for the current 



/ = dX deTr 



5i(e,A)/?o(e,A)L_i(e + ey) + 



P_i(e, X)Ko{e, A)Li(e - eV) + Xo(e)io(e) 



(36) 



which demonstrates that our numerical procedure con- 
sists of performing the double integral in e, A and em- 
ploying the matrix recurrence relations p2p and p3p at 
each step of the integration. The matrices depend on di- 
mensionless parameters 7, eV/A, T/A, x and a2.z (the 
parameter ai is excluded by ai — 1 — q;2 — as). 

Let us now present some results of our numerical anal- 
ysis of the problem in questioiii^S. To begin with, it is 
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FIG. 9: (Color online) Current-voltage characteristics for 
0^2 = aa = 1/3, 7 = 1 and T = 0.1 A case. The values of 
X are 0, 1.3, 2 and tt (top to bottom). 
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FIG. 10: (Color online) The phase-dependent current in the 
symmetric case (02 = 02 = 1/3) for 7 = 0.5, T = 0.1 A and 
eV = 2.3A, 2.5A, 2.8A, 3.2A (bottom to top). 



satisfactory to observe that several important features 
established in the lowest order in tunneling survive ~ 
though possibly with significant quantitative modifica- 
tions - also within our non-perturbative analysis which 
includes the effects of MAR. E.g., in Fig. [7] we display 
the I-V curves evaluated at different values of the phase x 
for the value Gi about 4.5 times smaller than G2 and G3. 
All our results are normalized to the normal resistance 
of the structure Rn defined from Eq. (l27l) as 

^ = (as + a3)(l - a2 - as) (Gi + G2 + G3) . (37) 

Comparing the results presented in Fig. [7] with those 
obtained perturbatively, cf., e.g., Fig. HJ in both cases we 
observe peculiar gap-like features at eV = A -I- Ag in the 
I-V curves at voltages A < eV < 2A. These features are 
also qualitatively consistent with the resultsi^ obtained 
for a model of Y-shaped Andreev interferometer with a 
normal part being substituted by a single mode quantum 
wire. For instance, this simple model allows to predict 
current peaks at voltages eV = eA{x) and eV — A + 
where eV = ^Aix) is the phase-dependent energy 
of the Andreev subgap bound state. These peaks have 
the same physical origin as those found here at eV = 
Ag(x) (cf. the curves displayed in Fig. [7]at eV < A) and 
eV = A + Ag{x), except in our case the current peaks are 
modified both due to specific transmission distribution in 
our structure and due to the influence of the parameter 
7- 

A typical current-phase dependence is depicted in Fig. 
|S1 This curve can be qualitatively compared to one eval- 
uated perturbatively and presented in Fig. [5] at voltages 
A < eV < 2A. Though these curves differ quantita- 
tively, their qualitative behavior remains somewhat sim- 
ilar. Namely, in both cases we observe an increase of the 
current with x a-t smaller phase values followed by its 
decrease at larger values of x- 

On the other hand, there also exist significant differ- 
ences between the results obtained in MAR and tun- 



neling regimes. Perhaps the most essential one is the 
presence of non-zero subgap current observed within the 
non-perturbative regime even in the limit of low tem- 
peratures. Obviously this feature is lacking within the 
lowest order perturbation theory in tunneling. The sub- 
gap current is well pronounced at Gi « G2 « G3, as 
it is demonstrated in Fig. [HI Similarly to the case of 
SNS junctiong^^-, there exist somewhat chaotic phase - 
dependent jumps of current in the subgap voltage region. 
At X = TT the current is strongly suppressed for subgap 
voltages eV < A. 

Our results demonstrate that MAR may essentially in- 
fluence the phase dependence of the current. In this re- 
spect it is instructive to compare Figs. [TUlandlSl While 
in the perturbative tunneling limit the current modula- 
tion decreases rapidly with increasing voltage (see Fig. 
[5]), there exist a clear voltage-independent modulation of 
the current in the non-perturbative MAR regime (Fig. 
[To)) . This modulation is due to the presence of a phase- 
dependent excess current with the amplitude ~ A/ [cRn) 
which is not captured within the lowest order perturba- 
tion theory in tunneling. 

Fig. [TT] illustrates the dependence of the current mod- 
ulation on the electron-hole dephasing parameter 7. This 
modulation clearly decreases with increasing 7 similarly 
to the minigap Ag, as it was discussed in the tunneling 
limit. It is worth pointing out that conductance asym- 
metry G2 7^ G3 also yields a decrease of the current mod- 
ulation, as it is demonstrated in Fig. 1121 

The phase-dependent current at lower voltages is dis- 
played in Fig. [131 Comparing these results to those in 
Fig. [6]we observe that the current takes much higher val- 
ues in the non-perturbative MAR regime. At the same 
time the current-phase dependencies turn out to be con- 
siderably smoother in this regime. Note that by tuning 
the voltage value one can reach the regime where the cur- 
rent depends monotonously on the phase Xi cf., e.g., the 
curve evaluated for eV = 0.5A. In this case the current 
modulation by the factor ^ 5 is observed. This bias volt- 
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FIG. 11: (Color online) The phase-dependent current in the 
symmetric case at eV = 2.5A, T = O.IA and 7 = 0.5,2, 10 
(top to bottom). 



0.30: 

0.25 : 




FIG. 13: (Color online) The phase-dependent current of the 
symmetric junction at T = O.IA and 7 = 1. The voltage 
values are 0.8A, 0.5A and 0.2A (top to bottom). 
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FIG. 12: (Color online) The phase-dependent current at eV = 
2.5A, T — O.IA and 7 = 0.5. The lower, middle and upper 
curves correspond respectively to Q2 = Q3 = 1/3, to Q2 — 
0.444, 03 = 0.222, and to Q2 = 0.555, Q3 = 0.111. 



age regime can be conveniently employed for magnetic 
flux measurements at lower voltages. This regime ap- 
pears advantageous as compared to, e.g., eV = 0.8 A and 
eV — 0.2A, since in the latter cases there exist extended 
flat regions with nearly ^-independent current values. 



V. CONCLUSIONS 

In this paper we have constructed a general theory 
of Andreev interferometers with three superconducting 
electrodes. Our analysis revealed a large variety of dif- 
ferent regimes and features which can be studied experi- 
mentally and used for performance optimization of these 
devices. 

In the case Gi ^ 6*2,3 the I-V curves can be evalu- 
ated perturbatively in the barrier transmissions. In this 
limiting case an important role is played by the Andreev 
minigap Ag(x) which may cause strong dependence of 
the current on the phase x not very large bias voltages. 



as shown, e.g. in Figs. [5]and|ni On the other hand, cur- 
rent modulation decreases rapidly as eV increases above 
2A, see Fig. El 

If all barrier conductances are comparable, Gi k, G2 ~ 
G3 , it is necessary to go beyond simple perturbation the- 
ory in tunneling and include the effects of MAR into con- 
sideration. Features related to the minigap A^ (x) persist 
also in this case, cf., e.g. Figs. |l]and[71 but the curves be- 
come much smoother and the current signal larger. Sig- 
nificant current-phase modulation can be achieved both 
in subgap (Fig. [131) ^'i^d overgap (Figs. fT0llT2l) voltage 
regimes, i.e. both these regimes can be used for suc- 
cessful operation of Andreev interferometers with three 
superconducting electrodes. In the regime of large volt- 
ages this modulation is due to phase-dependent excess 
current which is not captured within the perturbative in 
tunneling analysis. 

It is also important to add that both in perturbative 
and non-perturbative regimes large values of the current 
modulation can be achieved provided the parameter 7 
defined in Eq. (|10p remains sufficiently small. This mod- 
ulation decreases drastically for large values of 7 and, 
hence, such values should be avoided in Andreev inter- 
ferometers. 

Also we would like to make a remark concerning the 
effect of current noise. Although this effect deserves a 
separate analysis, it is clear already at this stage that in 
the MAR regime it would be desirable to avoid working 
in the limit of low voltages eV ^ A, since in this case 
one could expect dramatic increase of current noise which 
could compromise the operation of Andreev interferom- 
eters. Theoretical analysis of current noise correlator 

S{u;) = J dre''^'' {SI{t)6Iit - r)) , «(t) = /(i) - (/) 

in the case of diffusive SNS junctions in the zero fre- 
quency Iimiti2ii2^ reveals that at low enough voltages one 
has S/{I) oc 1/V, i.e. the signal-to-noise ratio should 
decrease with decreasing voltage. This effect is directly 
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related to MAR. Indeed, an effective charge transferred 
at voltages eV « 2A/n (with integer n) equals to q — ne. 
As the number of Andreev reflections n grows with de- 
creasing voltage, the charge q grows as g ^ (l-|-2A/|ey|) 
and, hence, the ratio S/ (/) cx q grows too. We also re- 
mark that, as in the case of Andreev interferometers with 



a normal electrode^i, one can also expect to observe noise 
modulation depending on the phase difference x- How- 
ever, at this stage we do not expect that this effect could 
alter our conclusion about low voltage regime being pos- 
sibly problematic for successful operation of Andreev in- 
terferometers with three superconducting electrodes. 
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